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ABSTRACT 

A priori  error  estimates  for  Galerkin  methods  for  numerical  approximation 
of  the  coupled  quasilinear  system  for  c = c(x,t)  and  p = p(x,t)  given  by 
V • [a (x, c) {Vp  - y(x,c)Vz}]  = 0 , 

V-  [b (x,c ,Vp) Vc]  - u(x,c,Vp)*Vc  = 4>  <x)  , 

for  x e t e (0 ,T] , and  appropriate  Neumann  boundary  and  initial  condi- 
tions are  considered.  Equations  of  this  type  arise  in  models  for  the  miscible 
displacement  of  one  incompressible  fluid  by  another  in  a porous  medium. 
Estimates  for  both  continuous  time  and  fully-discrete  time  Galerkin  methods 
are  presented. 
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SIGNIFICANCE  AND  EXPLANATION 


r 


l 

i 
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Equations  of  the  type  stated  in  the  Abstract  arise,  for  example,  in 
models  for  the  miscible  displacement  of  one  incompressible  fluid  by 
another  in  a i>orous  medium  as  in  chemical  flooding  of  oil  wells  to  recover 
more  of  the  oil.  The  variable  c - c(x,t)  corresjwnds  to  the  concentra- 
tion of  tl»e  fluid  used  for  the  flooding  process  while  the  variable 
P * p(x,t)  corresponds  to  the  pressure  at  the  location  x in  the  medium 
at  time  t. 

This  pai>er  gives  error  analysis  and  error  estimates  for  certain 
numerical  methods  for  approximate  solution  of  the  coupled  quasilinear 
system  of  parabolic  partial  differential  equations  in  the  Abstract. 

Both  continuous-time  anti  fully-discrete  time  methods  are  presented  and 
ana  1 yzcd . 


L. 


The  i es|H'i)s  ibi  1 i t y for  t ho  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MKC,  and  not  with  the  authors  of  this  rejxirt  . 
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GALERKIN  METHODS  FOR  MISCIBLE  DISPLACEMENT  PROBLEMS  IN  POROUS  MEDIA 


Richard  E. 


Ewing1  and  Mary  Fanatt  Wheeler* 


1 . Introduction 

We  consider  the  numerical  approximation  by  Galerkin  methods  to  a problem  arising  in 
the  miscible  displacement  of  one  incompressible  fluid  by  another  in  a porous  medium. 

A set  of  equations  [6)  and  (9)  modeling  the  pressure  p(x,y,t)  and  the  concentration 


c(x,y,t)  is  given  by 

i_  fJjL.  fi£ 

3x  '■y(c)  l3x 


o-n  I fe*  (&-'<«>  If))  * £ If))  ' 


3x  x )y  y 


and 


(1.2) 


~ <«»<♦, D.c.Vp)  - uxc)  + ^ ( 4 ($ , D, c , Vp)  - uyc) 
+ ^ <»<♦, D.c.Vp)  |^>  + ^ (B(*,D.c,Vp>  |^> 


“ ♦(x.y)  ~ + q(x,y)c  (x,y,t)  , 

for  (x,y)  < fl,  t * J “ (0,T)  , where  the  reservoir  f)  is  a bounded  domain  with  boundary 

3fi.  Here  - kx<x,y)  and  ky  « ky(x,y)  are  the  permeabilities  in  the  x and  y 

directions  respectively,  u = u(c)  is  the  local  viscosity  of  the  fluid,  y • y(c)  is 

its  density,  and  —■  are  the  components  of  the  dip  angle,  and  q.  the  imposed 

»X  a y 

total  flow  rate,  is  a linear  combination  of  a finite  number  of  Dirac  measures,  (i.e. 
sources  and  sinks).  The  function  if  = <t>(x,y)  is  the  porosity  and  D * D(x,y)  is  the 
molecular  diffusion  level.  The  components  of  the  dispersion  coefficient,  a,  4 and  B, 
are  functions  of  porosity,  diffusion,  concentration,  and  the  gradient  of  pressure.  It 
is  not  uncommon  in  petroleum  reservoir  engineering  for  one  to  assume  that  B * 0 and 
a - u(D)  and  <5  * 4(D)  (see  (9)).  The  function  c is  equal  to  the  specified  concentra- 
tion at  injection  wells  (sources)  and  to  c(x,y,t)  for  production  wells  (sinks).  No  flow 
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ooiut  i t ion*  jii«*  axmuiuM  on  t hoMifttat  v of  fit*  tonotvoii  awA  an  Initial  »*v»»w«*ih  i .it  ion  m >)lvr»n 


r 

1 


throughout  thr  rraervotr.  For  a grnrral  dtaottaaion  of  thr  phyal.n  of  mtaolblr  dtapla.r- 
ntmt  problrma  tn  ipafi vo 1 1 mg  marring  lh«  t radar  la  mfrrrrd  to  I'raorman  (7), 

lit  thr  ualrrkln  I'lwalm p wp  aha II  ua«>  diffrmnt  aubaparra  of  ll'(|!)  to  appro* t mat r 
pmaauir  amt  cenoant  rat  Ion.  w.<  not  a that.  If  r ( X , y . t ) 1,  than  (1.1)  amt  (l.J)  am 

tdrntlral.  In  thta  rasa  w.  am  fotvrd  to  uar  t hr  aamr  approx Iwat  lug  aubapaora  for 
pmaauir  an.t  conotmt  rat  ton.  Wr  air  thua  l r>1  to  mplaoltig  (l.J)  by  tho  non-divn  grn.-n  form 


whloh  ta  obtalnrd  by  multiplying  (1,1)  by  o amt  aubt  raot  tng  tltr  rrault  fmat  (l.J). 

Sat  tail,  l'ttv-r,  ami  lHt|X'itt  |s)  prrarntrd.  without  analyaia,  a finitr  rlmtrnt  mrth.st  rot 
solving  (l.l)-(l.tl  wttli  S 0.  for  alnt|'l  toity,  wr  shall  t mat  thr  oottpl  rd  guaa t I t nrar 
ayatrat  with  dri'rndrnt  vat  tab  Ira  o ■ nx.t)  an,t  p - p(x,t'  glvrn  by 


(1.4) 

V- 1 a ( x . o ) (Vp  - 

> (x,r)V*l|  - 0 . 

(l.S) 

V- |b(x,o,Vp)Vo] 

| - u(x,o,Vp) -Vo  - *(x)  , 

for  x » il,  t * J 

(0,T|,  whrm  t< 

la  a boumtrd  domain  in  , 

d - with  boundary 

,>t)  amt  u(x.r.Vp)  - •a(x.o)  (Vp  • t(x. 

,0)7*1  is  a vrotot  in 

Thr  analyaia 

for  (1.4)-(l.M  raatly  rxtrinln  to  thr  ,-a»p  (l.l)  amt  (l.t)  wham  mtxrd  x - y dm  tva- 
ttvra  am  ptrarnt.  t*»  aaauntr  that  thr  following  boundary  amt  initial  romttttona  hol.ti 


(l.h) 

- ><*.c)  JjJ  - 0. 

X 9 

.it), 

t r 

(1.7) 

b(x,c,Vp)  • 0 . 

X 9 

JO, 

t r 

(1  .It) 

P(X.O)  - l'0(x)  , 

X 9 

>;  , 

whrtr  '**  ta  thr  normal  dm  tvat  tvr  of  f on  thr  boundary  of  t!.  Wr  notr  that 

!)\> 

(1.4)-(l.M)  will  drftnr  p(x,t)  only  to  within  an  arbitrary  oonatant.  Wr  shall 


I 


normalise  p by  the  condition  that 


(1.9)  y^j-  / p(x.t)dx  - 1,  t « J , 

where  |i)|  la  the  measure  of  il.  Wo  notp  that  for  (1.4)  and  (1.5),  we  have  chosen 
q 0 in  (1.1)  and  (1.3).  If  q is  assumed  to  be  smooth,  which  corresponds  to  smoothly 
distributed  sources  and  sinks,  the  analysis  follows  with  few  chanqes. 

The  paper  consists  of  three  additional  sections.  In  Section  2,  the  continuous 
time  and  discrete  t isie  Oalerkin  procedures  are  formulated.  Assumptions  on  the  coefficients 
corresponding  to  the  physical  problem  and  smoothness  assumptions  on  the  solution  are 
given.  Methods  are  presented  which  allow  the  approximations  for  concentration  and 

pressure  to  lie  in  different  subspaces  of  II1  (kl)  . In  Section  3,  a-priorl  error  estimates 

2 

for  the  continuous  tins-  approximations  are  qiven.  Optimal  L rates  of  convergence  are 
established  for  the  case  where  the  dispersion  coefficient  depends  only  upon  molecular 
diffusion  using  subspaces  of  the  same  order  of  approximation.  In  the  case  that  the 

2 

dispersion  coefficient  also  depends  upon  the  Darcy  velocity,  wo  obtain  an  optimal  L 
rate  using  a subspace  of  one  higher  order  for  the  pressure  than  for  the  concentration. 

In  Section  4,  these  results  are  extended  to  the  discrete  time  case. 


I 


2.  The  Fmif  El— nt  Hsthods  and  • Precis  of  ths  Result! 

Let  (u,v)  • / uv  dx  and  ||u||^  - (u,u) . Let  W*(Q)  be  the  Sobolev  space  on  l) 

a * 


with  nor® 


y ll^ir 


% U-Ilk  • 

with  the  usual  modification  for  s » ■».  When  a - 2,  let  ||*||  • ||*||  II  ♦ II  w * lf 

i "a 

Vf  - (f,.f,)  write  ||  Vf  ||  in  place  of  (||f  ||S  ♦ ||f,||S.)*.  Also  H*  (30)  will 

W WK  Wk 

* s s 

denote  the  usual  Sobolev  apace  on  3(1 . 

Let  (.V^}  be  a family  of  finite-dimensional  subspaces  of  H1  (Si)  with  the  following 


property : 


For  p » 2 or  p * there  exist  an  integer  r » 2 and  a constant 
such  that,  for  1 < q < r and  * r w|’((})| 

inf  < II  ♦ - X II  0 * h ||  * - x II  . > I K0 II  ♦ II  hkl  . 

. - Li  ..A  , i'  ...g 


Similarly,  we  define  a family  of  finite-dimensional  subspaces  of  H til)  called  t.Vj  > 
which  satisfies  the  same  property  as  with  r replaced  by  s.  We  also  assume 

that  the  families  and  (A ’ ) satisfy  the  following  so-called  "inverse  hypotheses": 


< KQh  ||*||  - K0h 


b)  V* 


i _ lXph'l||v*| 

L (il) 


Restrict  il  as  follows  (with  S denoting  the  collection  of  restrictions) : 

2 

1)  rt  is  H -nKfular,  i.e.»  if 

-Av  ♦ 0v  • C#  x r . 0 - 0 or  1 , 


and  (C,l)  ♦ / ndjt  • 0 if  0 • 0 » 


IML  - Kdixlkll  ♦ llnll  >. 

H ‘(All) 

-)  ill  is  Llpschitc. 


4- 


b ; b(x,c,Vp)  5 b(x,c,p  ,p  ) 

x y 


u S u(x,c,Vp)  *=  (u^  (XiCi Pjj)  * u 2 ( x , c , Py ) ) . 

For  sow  t > 0 restrict  the  variable  q^  to  lie  between 

-e  £ 1 i + e . 

Assume  the  following  regularity  for  a,  y,  b,  u,  and 

(Q) : 1)  There  exist  constants  at,  b^,  and  such  that 

a)  0 * a.  i a(x,qi)  < Kx  , 

b)  | YJx.q^)  | i Kx  , 

c)  0 < * ♦ (x)  < K , 

(2.3) 

d)  | Vz  (x)  | < K , 

e)  0 < b,  <_  b(x,q,q2,q3)  , q'q2'q3  * K ' 

f ) lu^x.q  ,q  ) | < K (1  + |q  |>  , i « l,2,q4  « » » 

g)  |b(x,c,Vp) | < Kj  . 

2)  There  exists  a constant  M such  that  for  (c,p)  the  solution  to  (1.4)— (1.8) 

2 

and  for  (q2,q3>  t It  , i = 1,2, 

1 3a  3y  |a2b  |3u.  I 

|aic  (x,qi)  + <X>V  + — 2 (x'c'vp>  + j~  (x.q1.vp) 


dll . 

3 ui 

A *• 

3 u4 

3 u. 

dq2  (X'VV 

3x.3q,(x'c) 

X 2 

+ . (x,c,Vp) 

3c 

+ 

3c3q2  (X'C)| 

+ 1^7  <X'VW|  + ^X,q1,q2’q3)  | 1^7  <X'Wq3)|  • 

The  assumptions  made  on  the  coefficients  in  (Q)  are  physically  motivated.  In  the 

general  miscible  displacement  formulation  (1.1) -(1.3) 

2 2 


a = $ |d 


u 

X 

u 

1 r.  V 

lulll 

+ a 

fc  lllulll 

U U ^ 

* c»  ) 

_XJL1 

t' 

lllulll  > 

t - ♦ 


3 3 

I u u . 

Jo  ♦ a — ¥ — ♦ a — 5 — 1 , 

‘ III.. Ill  1 III. .Ill  I 


lull!  * |||u 

“here  HI  u HI  - J\i‘  » u*  an. 

x y 

cxn  shv'w  th.it  ci,  8,  aivi  5 satisfy  the  conditions  assumed  for  b,  that  K /u(c) 

x 

* /w(c)  satisfy  those  for  a,  and  that  u and  u satisfy  assumptions  made 

* y 


x -y  — 1 and  ai©  nonncgativ©  constants  (see  (6)  and  |9|)*  On© 


on  t h©  u . 

l 


l.©t 


♦<*.t) 


X"tq<a,b) 


1 < q < «• 


L ( <a.b) ,X> 

l*©t  (p»c)»  the  solution  of  (1 . 4)  - (1 . 8) , satisfy  the  following  regularity  assumptions: 


(R)  : 


(2.5) 


*»  IWI.  r*  £L  ♦ IIP  II  ..  siKi. 

L (JlH  ) 11  nL*(J|Hr)  L 2 


b) 


» ’«*  ' it 

t.  (.liH  1 " 


L (J;W  ) 


» c 

;«• 


L (J»L  ) 


• 1 

l.  <>t.H  ) 


I K;  < 


for  some  t > 0. 


The  analysis  proceeds,  following  Wheeler  1101,  using  a pair  of  auxiliary  elliptic 


problems . 

Let 

P«  ,vh 

<.\6) 

(a(  • .c( 

for  each 

t f J, 

where 

(J.7) 

1 

and  where 

<p.c) 

is  t he 

1 f 

iTT  1 pls.tldx  • 1,  for  each  t t J , 

M a 


ing  result. 

Leapta  3.1.  There  exists  K ^ » K , (C,  . K^. It  1 such  that 

Up  - pH  ♦ h||vip-  p>  || « K,h*||p||B  . 

Let  V > 3 be  chosen  sufficiently  large  that  the  bilinear  form 
- (b (c , Vp) V, . V\ 1 * (u (e , Vp) * V, , y ) ♦ \(,,y) 
is  coercive  over  ft' 01).  Let  c»  :f.  In'  the  elliptic  projection  of 


into  V , .let  inod 


- -(♦  . w)  ♦ 1 (c.w) , 


Vh  * 


t.'.gt 


l»(c,wi 


Plc.rfl 


w t 


? * ' h * * J*  Then.  a*  in  (10)  aiKt  (3),  w«»  ,411  obtain  the*  following  lumma, 

» Ttu*i*  exist*  K.  - K.(fi,b..V,K  ,K.,K  ,M)  suv'h  that 

4 4 ■ 0 1. 

I > (e  - c) 


*1 


L*  (J  iL  1 


L‘  (dll.*) 


♦ htlc  - c ||  , 

•>  I * IV.  U* 


L <J»H  > 


U" 


iVr{||c!!, 

1 L <JlH  > 


}• 


t.  (JiH  ) 


b,  ||^  - c||  ^ , < K.hr||c||  . . 

L (JlL  ) t.  (J;Hr) 

Ah  same*  that  there  exi»t»  a constant  *uch  that 

‘»*n>  ||vp||  _ _ ♦ II vc ||  „ . i «c,  . 

L (J.l  ) L <J»L  ) 5 

l-'  and  UO)  for  some  nu((  li'innt  condi  t ion*  fot  these  assumption*. 

We  tust  consider  continuous  t ime  approximations  of  ip,o).  Oenot e tin-  approxima 

t ton  of  p by  P 1 1 0 . T | • V and  t ho  approx imat  ion  of  o bv  p t 10,  Tl  » V , 

n h 

khoif  ti  v'l  is  defined  b\  t he  telations:  two  suppress  the  del's  tide  nos  ot  the  ooetffl 
. lent*  on  x' 


l-Wj  (a(C>Vr,Vv'  - (a  (0'>  1C)  V*,Vv)  , v » .V,  . 

h 

and 

l^.l"  (b  (0,  VPI  Vc,  Vw'  * (u  (0,  Vp)  • Vo.w)  * if  V • w)  * 0,  w t .V.  , 

at  h 

with 

1 C(x,0)  * o(x,0)  , 

w-ete  o t*  the  elliptic  project  ion  of  c defined  in  u'.'*'.  In  Sect  ion  I,  we  shall 
obtain  a priori  estimates  for  o - 0 and  for  V(p  - PI. 

Next,  we  define  a fully  disetete.  0(  ,\t ' -cot  rect  method  foi  approximat  ing  ip,c) 
to  ,si  on  backward  differencing  in  time,  l.et  At  ' 0,  N * T At  « J.- , and  t1'  - o.U . 
o < K.  Also,  let  j”  * i"ix'  * *(x.tnt  and  d(»n  - p,1’’1  - »■")  At,  Penote  the 
approximation  of  p by  w U'  - t'  . t ' tN  - r'  * \\  and  the  approx  imat  ion  of  o 


by  * « U ,t  1 * .V  . Assuming  that  wn  and  Z*'  art?  Known,  we  determine 

n 

*"  * and  zn  * as  follows: 

12. l*dt  t".t)  ♦ tbt*n.V^)7zntl.V*)  ♦ (u(Zn,\V')  •VZI'+i,*)  - 0.  »(  ,Vh  . 

And 

(2.1o)  (a(Zn  l)VWn  l.Vy)  « (a(zn  1 ) > i Zn* 1 ) V’r , Vy ) , y * N , 

h 

wtioie  ) • 0 or  I.  V>  note  that  the  coefficient  matrix  arising  from  the  algebraic 
system  (2.15)  with  J * 0 is  symmetric.  However,  in  many  problems  the  transport  term 
is  large  and  it  may  be  numerically  advantageous  to  use  (2.15)  with  j - 1 even  though 
the  coefficient  matrix  is  no  longer  symmetric.  Both  discretisations  are  0(At)  in 
time.  In  (2.1«>)  no  explicit  time  discretisation  error  is  made.  A priori  error 
estimates  will  be  presented  for  c"  - 2°  and  V(Pn  - wn)  in  Section  4. 

let  fn' 1 * (f1'*1  * tn)/2.  If  (2.15)  is  replaced  by 

1 1 
* u*  ~ 

(2.17)  <*d  i ,*)  ♦ (b(EZn,VEw")VZ  * , V*)  ♦ (u (Ez" , Vew")  • VZ  * , *)  - 0.  K>  t .V  , 

z h 

. -n  3 *n  1 *n-l 

***•*•  E*.  » - *.  - - Z , win  analysis  similar  to  that  given  in  Section  *4  (see  18], 

1 10],  12)  and  IS))  will  yield  a time  discretization  error  of  size  0((At)  j. 

When  estimating  the  actual  physical  quantities,  we  know  that  the  pressure  is  much 
smoother  in  time  than  the  concentration.  Tims,  in  practice,  one  should  not  use  (2.16) 
to  redetermine  VW  ^ at  each  time  step.  Instead  one  value  of  Vwn>*  should  be  used 
in  the  ^'efficients  of  (2. IS)  for  several  time  steps  before  a new  pressure  is  approximated, 
one  way  to  view  this  procedure  is  to  have  tvA>  different  time  steps  for  the  different 
equation*.  An  analysis  of  this  method  will  follow  in  a more  general  paper. 


^ 

1 

1 

j.  A I t ion  Error  k -t  mates  for  the  Continuous  Time  Approximations 

In  this  section,  we  develop  a priori  bounds  for  the  errors  (c  - Ci  and  V{p  - P) . 

We  shall  see  that  if  b depends  on  Vp  then  for  the  best  rates  of  convergence  one 
should  use  subspaces  of  piecewise  polynomials  of  different  degrees  for  {A/^}  and  ) . 

We  first  obtain  an  estimate  for  P - p where  P is  defined  by  (2 . 12) - (2 .14)  and  p 
is  defined  by  (2.6)-(2.7). 

We  assume  that  for  some  £ > 0,  c c,  (see  the  statement  preceding  Q)  , 

(3.0)  -t  < C < 1 + E . 

We  shall  prove  that 

|| c - C ||  - 0(hr  + h8"1)  , 

which  yields 

(3.1)  || c - C ||  _ « K(hr  1 ♦ hS  2)  . 

L 


Since  0 c 1 , one  can  deduce  from  (3.1)  with  the  restriction  r 2 and  s ^ 3 
that  (3.0)  holds  for  small  h. 

Lensna  3.1.  There  exists  a positive  constant  (a4 , , K^.M)  such  that 


|V(P  - p)  II  < k6I|c  - c | 


t f J . 

Proof.  Subtract  (2.6)  from  (2.12)  and  use  (1.4)  and  (1.6)  to  obtain 


(3.2)  (a(C)V(P  - p)  ,Vv)  - (la(c)  - a(C)  ]Vp,Vv)  ♦ ( |a(C)y(C)  - a(c)y(c)  )Vz,Vv)  , v ( *h  . 

If  v » P - p € N.  t then 
h 

aJ|V(P-  p)  ||  2 < | (a  (C)  V (P  - p),V(r  - p) ) | < ||c  - c||  ||  V (P  - p)  || 


(3.3) 


iy  ||V(P-  p)  II2  ♦ 


c-  c| 


t f J , 


where  K depends  upon  uniform  bounds  for  |^-|  , ilii.  |Vz|,  and  || Vp ||  ^ ^ . The 

Sc  dc  L (J;L  ) 

result  of  the  lemma  follows  directly  from  (3.3). 

We  next  compare  C and  c,  defined  in  (2. 12) - (2 . 14)  and  (2.9),  respectively.  We 

shall  make  an  additional  assumption  on  P defined  in  (2 . 12) - (2 .14) . Assume  that  there 

* 

exists  a positive  constant  K such  that 


(3.4) 


iVPH  . - i K 

L (JtL  ) 


Without  loss  of  generality,  assume  K > 2KS. 
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Theorem  3.2.  There  exists  Xg  - Kfl <* >b, >♦* #MiKi , i < 6;K  ) such  that 


0.5)  I|C  - e||  „ , + ||V(C  - c)  ||  < K (hr  + h8_1)  . 

L (J)L  ) L (J|L  ) 

If  s and  r from  Lemmas  2.1  and  2.2  satisfy  s > 3 and  r > 2 and  h is  taken 
sufficiently  small,  then 

‘3.6)  IMI  . . < 2K5  < K* 

L (JiL  ) 

* 

and  a Kfi  above  can  be  chosen  which  is  independent  of  K . 

Proof.  Subtract  (2.9)  from  (2.13)  to  obtain 

(b(C,VP)V(C  - c)  ,Vw)  + - ( (b(c,  Vp)  - b(C,VP)]V^,Vw) 

(3.7)  +(*(ft‘ft)'W  - (u(C,VP) -V(C  - c)  ,w) 

+ ( [u(c,Vp)  - u(C,VP))-Vc,w)  + X (c  - c,w),  we  M.  . 

h 

Next,  letting  w » C - c e M , integrating  the  left  side  of  (3.7)  from  0 to  t and 

n 

using  (2.14),  we  obtain 

t . t 

/ (♦  T-  (C  - 8),(C  - c)  )dt  + / (b(C.VP)V(C  - c)  ,7  (C  - c)  )dt 

0 3T  0 

(3.8) 

>7  ||(C  - c)  (t)  II 2 +b , J II  V(C  - c)  (x)  II2  dl  . 

0 

We  next  use  (2.8),  Lemma  2.2,  Lemma  3.1  and  Holder's  inequality  to  see  that 
t 

|/  ((b(c.Vp)  - b(C,VP))V8,V(C  - c) )dt  | 

0 

t 

^ K l|Ve||  „ „ / ( ||  c — C ||  + ||  V(P  - P)||)||V(C  - c)||di 

L (JiL  ) 0 

t 

1 Kln  / (||e-  a||  + ||e-c||  + ||  V (p  - P)  || ) ||  V (C  - c)  ||dT 


* Of  2 

<—  / ||V(C-  C)  (T)  II  dT  + K / ||(C  - C)  (T)  II  dt 


+ K / (||  (c  - C)  (T)  |r  + II  V(p  - p)  (T)  II  >dT 

0 

< -T-  / II  V(C  - c)  ( T ) II  2dT  + K / II  (C  - c)  (1 ) ||2dl  + K (h‘  r * h‘  S 2)  . 
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MffOMEWMNMMEW 


Here  depends  upon  M,  Kj  , K^,  , and  bt,  while  depends  upon  M,  K(  , K^,  K( 

K^,  K^,  Kfe  and  b#.  We  can  similarly  obtain  the  bound 
t t 

(3.10)  |/  (lu(c.Vp)  - u(C.VP)  | -Vc.C- ^)dT|  < K / ||  (C- c)  (t)  ||2dt  + K (h‘Jr  ♦ h2‘_2)  . 

0 0 1 


Next,  from  (2.3),  we  note  that 


(3.11) 


t 

1/  (u(C.VP)  -V(C  - c)  ,C  - £>dt| 

0 

- Klb(  II VPH  ->  .»  ♦ 1,3  / II  (C  - ^)(T)||2dT4  ~ ! ||V(C  - £)<r)||2dT  . 

1.  <J|L  ) 0 H 0 


Then  using  l.emma  2.2,  we  obtain 


(3.12)  |/  (♦(ff- |f)  - C- ojdr  * X / (e-c,C-c)dT|  < K1?  / ||  <C  - £>  (t)  ||  2dr  + K1?h2r  . 

Combining  like  terms  in  the  estimates  0.8)  - (3.12)  and  the  assumption  (3.4),  we  obtain, 
from  (3.7), 


Y ||  (<’  - c)  (t)  ||2  + y / || V(C  - c)  (x)  ||2dt 
(3.13)  0 

l*ls  / II  (C  - c)  (t)  ||2dt  3 K19(h2r  3 h28'2)  , 

* 

whete  depends,  in  particular,  upon  K . Then  applying  Gronwall's  lemma  to  (3.13), 

we  obtain 


( 3.14) 


lie 


‘H  . a 

1.  (Jit  ) 


||  V (C  - c)  ||  < K (hr  3 h8'1)  , 

L (J|L2) 


where  we  note  that 

<3-*5>  K20  < Knexp(K22K*2)  . 


To  complete  our  argument,  we  must  allow  that  for  h sufficiently  small  (3.e)  holds. 
We  use  lemma  1.1,  (2.2),  (2.11),  (3.14),  and  (3. IS)  to  see  that 
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M „ . < II v tr - p)  ||  _ a,  ♦ ||vp||  ^ . 

L (J|L  ) L (J|L  ) L (JiL  ) 

iK  ♦ K h'1||V(P-  p)||  „ 

I.  (JiL  > 


(3. lt>) 


i Ks  * V’~1k6<I!c'  «ll 


♦ c - c 


L (Jit  ) 


» 2 

L (J|  L ; 


- K5  + K0K6K2i'**r(K22K*2'<*'r'1  + **S'2>  ♦ KoK2K4V'r'1 


Thou  clearly  if  l\  is  taken  sufficiently  small,  and  If  r > 2,  and  s - Ji  then 


(3.17) 


II  VP  || 


I.  (J It  > 


<_  < K . 


Wo  see  that,  from  (3.5),  in  order  to  got  the  best  rates  of  convergence,  one  should 

use  a space  of  piecewise  polynomials  of  degree  one  qreater  for  N.  that  for  .V.  . For 

h h 

example,  the  use  of  piecewise  linear  polynomials  (r  - 2)  for  and  piecewise 

quadratic  polynomials  (r  - 3)  for  . would  yield  a rate  of  convergence  of  the 
2 

form  0(h  ) for  C - c. 

We  next  combine  the  results  of  Lemma  3.1  and  Theorem  3.2  with  Lemmas  2.1  and  2.2 
and  the  triangle  inequality  to  obtain  the  following  result. 

Theorem  3.3.  There  exists  a constant.  K,(  such  that 

(3.  IB)  h ||  V (P  - p)||  „ > ||C-  C ||  iv  ♦ h||  V (C  - c)  ||  < K (hr  ♦ h*'1)  . 

L (JlL  ) L (JiL  ) L (J|L  ) 

Some  mathematical  models  for  miscible  displacements  which  are  currently  being  used 

by  oil  companies  make  the  assumption  that  the  coefficient,  b in  (1.2)  does  not  depend 

uivn  Vp,  but  only  upon  x and  c.  Making  this  assumption,  we  obtain  optimal  order 

2 

convergence  tales  in  the  I,  norms  using  test  spaces  composed  of  piecewise  polynomials 
of  equal  degrees.  These  results  involve  deriving  a sharpet  estimate  than  obtained 
in  (3.10). 

Theorem  1.4.  Assume  h - b(x,c)  in  (1.2)  is  independent  of  Vp . There  exists  a 

« 

constant  K,.  K , . ( V ,b.  ,4  . ,K  K ,i  - e;K  ,M)  sneh  that 

24  24  * * 0 i 

„ , ♦|l°'v||  v , * h|)  V(C  - c)  ||  , - K <hr  » h*  • h*'”"’' 

I (J;  I.*)  1.  (.1:1  ) l‘(Jil.  ) 


-1. 


I ».  I")  h|| V (P  - p)  || 


wm* 


If  either  r > 3 or  s > ] for  r ■ 2 and  h is  taken  sufficiently  small,  then 


(3.20) 


VT 


- „ i «5 1 K 

I.  (J|L  > 


and  a K_.  can  be  chosen  which  is  independent  of  K . 

24 

Proof . Most  of  the  estimates  will  follow  as  in  the  proof  of  Theorem  3.2.  Let 
C - C - c.  The  estimate  (3.9)  can  be  replaced  by 
t t 

1/  ((b(c)  - b(C)lVc,Vc)dt|  < KJ5  / ( ||c -i||*  ||C||)  ||vc||ar 


(3.21) 


lY  ! I|vc(t)||2dt  * k25  / II  C(T>  ||  2dT  ♦ K26h2r  . 


bM  t 

/ 

0 ‘■'0 

In  order  to  obtain  a sharper  estimate  for  (3.10),  we  rewrite 

<[u(c,Vp)  - u (C,  VP)  ) • Vc , c)  - ((u(c,Vp)  - u(c,Vp)  1 *Vc>0 

(3.22)  ♦ ( [u  (c,  Vp)  - u(c,Vp)J.V(c  - c),C) 

♦ Uu(c.Vp)  - u(C,VP) ) *Vc,C)  » T ♦ T2  ♦ T3  . 

3u  3u 

Prom  Section  1,  we  note  that  for  u - u(c.q),  - ~-i  (c)  is  independent  of  q,  i-  1,2. 

Integrating  the  first  term  on  the  right-hand  side  of  (3.22)  by  parts,  we  have 


f - 3 

r Jui 

3c  1 

3 r3u2 

3c  1] 

(p  ’ p'  *7 

- — (c) 

L*! 

3*!  "J 

* 3*2  I-9'); 

(C>  ^ CJJ 

(3.23) 


+ / (P  ~ P)C 
3(1 


3u  , 3u,  , 

3q~  (c)  3^  vi  + 3q7  (c)  aTT  V2 


'1  * 2 

Thus,  from  the  assumed  smoothness  of  u and  c,  we  have 

Tsu 


ds  . 


(3.24)  ITjI  < K27I|P-  Pi 


♦ 


/ (P  - p)  c 
31) 


3u . 


1.  . , 3c  ~~2  . . 3c 

^7  <c>  3^  V1  * 317  s 


vis 


We  now  estimate  the  above  boundary  term  by  considering  the  following  Neumann  problem: 

a)  -V-a(x, c)Vi|'  * i/'  ■ 0,  x r a , 

ait 


(3.25) 


b)  a (x,c)  “ - p 


x r 3fi 


r " (c)  £7  V1 + w, <c)  £7  v2] 


By  our  regularity  assumptions  (S)  and  (R)  , we  have  i|i  t H (0)  and 

(3-27)  ll*||,  £ K ||  T || 

2 II  (90) 


We  note  that  since  c * L (Jjh  (0))  and 


a(c)  is  uniformly  bounded, 


G1  t L (J;H1//2+<  (30))  and  we  have  by  Lemma  2,2  of  [1], 

Ill’ll  1/2  £ K II Gjll  „ / lie 

,,  (30)  L (JjH  (30))  H*'‘(30) 


Multiplying  (3. 25. a)  by  p - p and  integrating  by  parts,  we  see  that 

(a (x,c) V<i, V (p  - p))  t (ip,p  — p)  — / l'(p-p)ds  “ 0 . 

30 

Using  (2.6),  we  note  that 

(3.29)  / r (p  - p)ds  - (a  (x,c)  V(i|j  - i|i  ) ,V(p  - p) ) + (i|i,p  - p)  , i|i  e N.  . 

30  h 

Hence,  by  (2.1),  Lemma  2.1,  (3.27),  (3.28),  and  (3.29),  we  have 

1/  r(p-p)ds|  £ K2  (h||v<p-p)  ||  + II  p-p  II  ) ||||)||  £ K (h||v(p-p)  ||  t ||p-p||  > ||  T || 


1 1/2 
H (30) 


£ K30{hllV(p"P)  II  + IIp-pIIHIcHj  <~  ||t||j  + K^h23  . 


Combining  (3.24),  (3.30),  and  Lemma  2.1,  we  obtain 


lTll  - T l|c|lJ  + K32ll^l|2  + K33h2 


Using  Lemmas  2.2  and  3.1  we  note  that  the  third  term  on  the  right  of  (3.22)  can  be 
bounded  as  follows 

(3.32)  |T  | £ K ||  Vc||  „ „ (||V(p-  P)  ||  + K*  ||c  - c||)||c||  < K ||  C ||  ^ + K h2r  . 

L (JjL  ) 


In  ordar  to  bound  the  second  tern  on  the  right  eide  of  (3.22),  we  use  (2.2)  and 
2.1  and  2.2.  We  obtain 

" ' “ Id 


(3.33) 


It2I  1 K l|V(p-p)  II  „ II V (c  - c) 

L (J)L  ) 


£ K39h-V-Mld|  £ k40IUI!2  ♦ K41fc' 


L 

2r<2s-b 


£ Kjgh3”^ || V (c  - c)||  ||c| 


Usinq  the  estimates  (3.8),  (3.11),  (3.21),  (3.22),  (3.31),  (3.32),  and  (3.33)  in  (3.7), 
we  obtain 

♦ . „ ib.  t b„  t - t 

(3.34)  2 


U<t>H2+-J-/  ll’C(T)  ||2dT  £ ~ I II  C <T>  ||2dt  + K42  / ||c(T)||2dT 


♦ KaJh2r  ♦ h2S  ♦ h2r+28-6)  . 

4 J 


We  can  next  add  — — / ||c(t)||‘dT  to  both  sides  of  (3.34)  and  use  the  definition  of 

, 0 

the  H norm  to  hide  the  first  term  on  the  right  of  (3.34)  on  the  left  side  of  the 
resultinq  inequality.  Then  an  application  of  Gronwall's  lemma  yields 

(3.35)  lldl  „ , + ll»t| 


L ( J ; L ) 


I 2 2 — *44^  + ^ + hr+S_3}  • 

L (JjL  ) 


The  rest  of  the  theorem  follows  as  in  Theorem  3.2  and  Theorem  3.3  with  the  minor 
modification  of  usinq  (3.35)  instead  of  (3.14)  in  (3.16)  and  the  assumption  that 
r ♦ s > 5. 
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t 


4 . A Priori  Error  Ksti mates  for  the  Discrete  Time  Approx imat ions 

In  this  section,  we  develop  a priori  bounds  for  the  errors  c"  - z”  and  V(pM  - Wn 

at  discrete  time  levels  tn  where  z”  and  W°  are  defined  in  (2.15)- (2.16) . The 

results  obtained  will  be  similar  to  the  corresponding  estimates  for  the  continuous 

case  from  Section  3.  We  assume  that  S,  Q,  and  R and  the  restrictions  on  {3/^ } and 

( S , } of  Section  2 hold, 
h 

Theorem  4.1.  Let  Z°  be  determined  such  that 

*>  I|Z°  - cjl  i K45hr 

(4.0)  and 

b)  -«  < Z°  < 1 ♦ * . 


The 


i>re  exists  a constant  K « K (1  ,b.  ,M;  K , , i < 6)  and  constants  t > 0 and 

46  46  * * i — 0 


h^  > 0 such  that  if  At  <_  t and  h hQ  , 


N-l 


sup 


|Z  - c ||  ^ 4 h2  J ||  V (2  - c)n||2At  ♦ h2  sup  || V (W  - p) 


n»0 


(4.1) 


< K ( (At)2  4 h2r  4 h2(S"U  ) . 

— 46 


Proof . Assume  that  Z satisfies  the  induction  hypothesis 


(4.2) 


- 1 < 2 < 1 4 f 


for  n » 0,l,...,t  - 1.  Let  c"  ” Zn  — c"  and  n"  ■ Wn  - pn.  Subtracting  (2.6)  from 
(2.16)  we  have,  for  n * 0,1,..., N - 1, 


(4.  1) 


, „ . ,,  , n+1.  ,_n+l , , „~n4 1 „ , 

(a (Z  )Vn  ,Vy)  - ( (a (c  ) - a(Z  ))Vp  ,Vy) 

. ,,  ...ntl,  , ntl  , n4l,,„  _ , 

4 ( la  (2  ) Y (Z  ) - a (c  )y(c  )]V2,Vy), 


V * #h  • 


As  in  l emma  3.1,  we  let  y ■ * N , use  (4.2)  and  (tf)  to  obtain 


(4.4) 


l|Vn"+1||  < K47(||tn+l||  4 hr). 


n « 0, 1 .... ,N  - 1 . 


Next,  subtract  (2.9)  frivn  (2.15)  and  use  the  notation  from  Section  2 to  obtain,  for 
n - 0,1,..., n - 1,  and  .1  = 0 or  1, 
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(4 . S) 


X) 


(♦dttn,x)  ♦ (b(Zn,Vwn)Vcn+1.Vx)  - (t""1)  -d^"], 

- X ( (c  - £>n+\x>  ♦ (lb(cn+1,Vpn+1)  - b(Zn,Vwn)]Vcn+1,Vx) 
♦ (lu(cn+l,Vpn+l)  - u(zn,VWn)] -Vcn+j,x) 


- <u(zn,vwVvcn+j.x>. 


For  simplicity,  let 
(4.6) 


X * * 


Now  lettinq  x “ t”*1  * and  using  (4.6),  we  estimate  the  left  side  of  (4.5) 


by 


~ (*(C"+1  - tn).tn+1)  ♦ (b(Zn,VWn)Vcn+1,Vcn+1) 


(4.7) 

- alt  { II cI1+1  II*  - IUnll^>  ♦ bjvcn+1||2  • 

We  next  use  Lemmas  2.1  and  2.2,  (2.11),  and  (4.4)  to  see  that 
|([b(cn+1,Vpn+1)  - b(zn,vw") ivcn+1,vcn+1) | 

- | (lb(cn+1,Vpn+1)  - b(cn,Vpn+1)  + b(cn,Vpn+1)  - b(zn,Vpn+1) 


(4.8) 


♦ b(zn,Vpn+1>  - b(zn,vPn)  + b(zn,vPn)  - b(zn,vwn)ivan+1,vcn+1)| 

< K50(6t||dtcn||  ♦ ||(c-  c)n||  4 |kn||  > At||dtVpn||  4 ||V(p-  P)n|| 

♦ l!vnn||>  -||vtn+1|l 

i K51||cn||2  4.  ^ ||VCn+1||2  + K52(At)2{||dtcn||2  4 ||dtVpn||2)  4 K53{h2r  4 h2*8'1’}  . 

In  a similar  fashion,  we  obtain 

|<lu(cn+1.Vpn+1)  - u(Zn,Vwn)).Van+j,Cn+1)| 

(4-9>  I ~ Ill'll2  ♦ K54Hcn||2  + K„(At)2{||dtcn||2  4 ||dtVpn||2) 

4 K (h2r  4 h2'3'1')  . 

56 


Next,  from  <2.3  f) , we  note  that 


(4.10)  |-<U<Zn,VwVvtn+’,!.nn)|  i K57(||vWn||  ^ 4 K5e)2||rn+1||2  4 ~ ||VCn+j||2 


L (Si) 
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We  then  estimate  the  first  two  terms  on  the  right  side  of  (4.5)  using  l.emraa  2.2  as 
follows i 

l,.,3c  . n . . n , -n,  n41,  ...  - n+1  n+l% i 

* laT  (t  ' Jtc  dtc  ' dtc  ) - M(c  - c)  ,t  )| 


(4.11) 


where 


(4.12) 


< £ IUn+1ll2  ♦ V»(c-  *>n+lH2  ♦ ll«t<e-  ='n!t2'  * *60°2 

II  "1H2  + Khih2r  + Kh00n- 


tn+l  \2 

i |Nc.»)|k 

n "at*  11  J 


tn+1  ||  _ 2 

I At  / — ^ ( • ,s)  ||  ds 

tn  M3t 


by  Schwarz's  inequality.  We  note  that 
N-l 

(4.13)  At  l & <_  (A 
n-0  n 


N-l  tn+1 

U2 

2 it  ,2  , 

l i 

r§  (,-a) 

ds  < (At)"P— 

n-0  ^n 

'at 

Nat21 

L2(J|L*) 


We  next  multiply  each  of  the  estimates  (4.7)-(4.11)  by  At  and  sum  on  n,  n - 0,...,<  - 1 
to  obtain  from  (4.S) 

f-1 

l 

n»0 


I (♦.lie'll2  - ♦*||C°||2)  ♦ b.  I l!vCn+,||2At 


i T j (iun4,ii5  - iunu2>  > b.  Y iivcnnn2At 

n-0 


n-0 


(4.14) 


Kb2  T lkn||2At  ♦ ~ Y ||Vcn,1||2At  4 At  £ lie 

n-0  — " 


b^  t-1 

y 

n-0 


* II . f li  2 


f-1 


♦ Kc,  l <I|VW"||  ■*  K.  >‘||cn4l||2At  + K , ( (At) 2 ♦ h2r  ♦ h“  (*-X)  > , 


57 


OV 


n-0  I.  (i!) 

We  see  that  an  application  of  Gronwall's  Itmi  would  ooisplete  our  argument  if  we  did 
not  have  the  next  to  the  last,  term  on  the  right  of  (4.1s).  We  shall  use  an  induction 
argument  to  treat  this  term.  As  an  induction  hypothesis  assume  that  toi  ti  sufficiently 
smal  1 


(4.15) 
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l (l!) 


< 2KS>  for  n - 0,1 - 1 . 


Usinq  (2.2),  Lemma  2.2,  (2.11),  (4.0),  and  (4.41  as  in  (3.1e)  wo  see  that  (4.1r>)  clearly 
holds  for  t ■ 1.  Using  (4.15)  wo  see  that  if 

(4.16)  it  < ♦,(4K57{2Ks  ♦ K }2)-1  = t0  , 


then  wo  have 

(4.17)  |kV‘b.  V ||Vtn4 'll  *’ At  < Y ||  t*'  ||  * At  ♦ K ((At)2  t h2r  * h2*8'1’)  . 
n-0  n«0 

Then  an  application  of  dronwall's  lemma  gives  us  an  estimate  which,  with  (2.2),  l.emma  2.2 
(2.11),  and  (4.4)  shows  as  in  (3.16),  that  the  induction  hypotheses  (4.2)  and  (4.15) 
holds  for  n * t.  finally  an  application  of  the  triangle  inequality  coupled  with  hemmas 
2.1  and  2.2  and  (4.4)  yields  the  desired  result. 

.lust  as  in  Section  3,  if  we  make  the  assumption  that  b ■ b(x,c)  in  (1.2)  is 

2 

independent  of  Vp,  we  can  obtain  optimal  order  L rates  of  convergence.  The  proof 
will  follow  by  combining  the  techniques  of  Theorem  3.4  with  those  of  Theorem  4.1.  We 
obtain  the  following  result. 

Theorem  4.2.  Let  the  assumptions  of  Theorem  4.1  hold.  Assume  b - b(x,c)  in 

(1.2)  is  independent  of  Vp.  There  exist  constants  K._,  t._,  and  h,  such  that,  it 

66  0 0 

At  ^ tQ  and  h <_  hQ, 


sup  ||  Z-  c ||  * + h T 1 1 V ( 7.  - c)n||‘At  + h*  sup||V(W-  p)  | 
. n n»0  n 


< K ((At)2  ♦ h2r  ♦ h29  ♦ h2r+2»-6)  . 
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As  we  mentioned  at  the  end  of  Section  2,  if  we  replace  (2.15)  by  (2.17)  and  extrap- 

2 4 

olate  the  coefficients,  we  can  replace  the  (At)*  in  (4.17)  and  (4.1H)  by  (At)  and 

2 

get  discretisation  errors  in  time  which  are  CM  (At)  ).  In  order  to  do  this  we  must 
determine  a starting  procedure  to  obtain  the  approx imation  at  time  t'  - At  since 
two  levels  must  be  known  to  determine  the  next  level  with  this  method.  A predictor- 
corrector  version  of  (2.16)  will  suffice  as  a starting  procedure  of  sufficient  accuracy. 
Since  the  proof  techniques  are  similat  to  those  presented  in  Theorem  4.1  and  are 


-in- 


»imUt  to  those  |'i  evented  tu  111.  Ul.  IS).  |H|  and  |10|.  we  shall  not  present  thru.  Imr. 
w»-  luquiie  slightly  »oi  e smoot hues*  on  <-  and  p for  these  result*. 


rtiroirm  4.  t.  If  wo  the  .tJJrJ  amoot hnr »h  on  e aik)  i'  that 


* Ob 


l.  (JiL  ) 


3t  " 2 2 

L (JiL  ) 


0 X 

«*ik1  if  .tiki  2 cur  Jetrrwu\<xt  to  »Atl»fY 


*'  II  *°  - e0l|  ♦ |l*1  - c(tX)||  a/  . 


b)  -»  » Z < 1 ♦ r , n - 0,1  , 

thou  there  mi»i  constant*  k , t , and  h.  such  that,  if  At  « t and  h < h , 

0 0 0 0 

> 2 N;l  „ 

"UP  II 2 - c||  ♦ h"  J ||V(K-  cl  II  At  ♦ h2  sup  ||  7(10-  p>  || 2 

( 4 . 1 ) t"  n“°  ,» 

' K {(At)4  + h"r  * h*"'* ) . 
o * 

finally,  in  the  ease  that  b * b(x,c) , we  can  use  the  technique*  of  Theorem  4..' 
to  til' tain  the  followin')  result  for  2 defined  in  ii.17). 

Tliev't.'i*  4.4.  Under  the  assumptions  of  Theorem  4.  J,  there  exist  constants  k t . 

70  0 

aixi  h^  such  that,  if  At  « t0  and  h ^h  . we  have 

N-l 

supllz  - clr  ♦ »'*  y ||  v (K  - c)n||*At  * h*  sup'lvtw-  p)||2 
, n n-o  n 


- K70UAt)4  * h*’r  * h*’B) 
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